Three-body resonant radiative capture reactions in astrophysics 
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We develop the formalism based on the S-matrix for 3-^3 scattering to derive the direct 
three-body resonant radiative capture reaction rate. Within this formalism the states, which decay 
only /predominantly directly into three-body continuum, should also be included in the capture rate 
calculations. Basing on the derivation, as well as on the modern experimental data and theoretical 
calculations concerning 17 Ne nucleus, we significantly update the reaction rate for lj O(2p,7) 17 Ne 
process in explosive environment. We also discuss possible implementations for the 18 Ne(2p,7) 20 Mg, 
38 Ca(2p,7) 40 Ti, and 4 He(na, 7 ) 9 Be reactions. 

PACS numbers: 21.45.+V - Few-body systems, 26.30.+k - Nucleosynthesis in novae, supernovae and other 
explosive environments, 25.40.Lw — Radiative capture, 25.40.Ny — Resonance reactions. 
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I. INTRODUCTION 

The reactions of the three-body radiative capture may 
play a considerable role in the rapid nuclear processes 
which take place in stellar media under the conditions of 
high temperature and density. The possibility to bridge 
the waiting points of the rp-process in the explosive hy- 
drogen burning by the 2p radiative capture reactions 
was discussed in Ref. The reactions 15 0(2p,7) 17 Ne, 
18 Ne(2p, 7 ) 20 Mg, and 38 Ca(2p,7) 40 Ti could be a more ef- 
ficient way to "utilize" 15 O, 18 Ne, and 38 Ca, than to wait 
for their (3 + decay (corresponding lifetimes are 122, 1.67, 
and 0.44 seconds). The 4 He(na,7) 9 Be reaction has been 
found to be important for building heavy elements in the 
explosions of supernovae [E Q ■ This reaction has been 
several times theoretically considered in the recent years 



The three-body radiative capture is a very improba- 
ble process. It can only be important if the sequence 
of the two-body radiative captures to the bound states 
is not possible. This happens if the bound intermediate 
system does not exist (along the driplines the contin- 
uum ground states are not uncommon). The three-body 
radiative capture can proceed sequentially via the inter- 
mediate resonances or directly from the three-body con- 
tinuum. The later is the inverse process to the "true" 
two-proton radioactivi ty JlOl- which studies are active 

now [II E3, EE 13 EE tlp^M EE El, Th e relations 
between sequential and direct mechanisms of two-proton 
decay are discussed in details in Refs. [l2llT8|. 

In the modern literature exists some misunderstand- 
ing about the role of the direct three-body capture in 
the theoretical calculations of the three-body radiative 
capture rates. In the most cases this misunderstanding 
does not lead to any significant problems. However, in 
some situations the difference is sufficiently large. In our 
opinion the origin of the misunderstanding is the follow- 
ing. The accurate formulae for the resonant three-body 



capture are known for a long time (see e.g. Ref. |22fl). 
But nowadays they does not seem to be always inter- 
preted completely correctly. The possible reason is that 
for nonresonant capture calculations the sequential cap- 
ture formalism is used 0, IE IE IE 113 • At some stage it 
has become considered as obtained in more general as- 
sumptions (see e.g. Ref. 23]) than the derivation of Ref. 
|22| , based on complete thermal equilibrium and detailed 
balance. 

In this paper we use formalism based on the S-matrix 
for 3 — » 3 scattering to derive the reaction rates for the 
three-body resonant radiative capture. In this approach 
the right way of using these formulae becomes evident. 
We find that the direct and sequential capture mecha- 
nisms are complementing each other. In the cases when 
the sequential process is prohibited energetically or sup- 
pressed dynamically the sequential formalism should un- 
derestimate the rate. Among the processes, where signif- 
icant differences with previous calculations can be found 
are reactions leading to 17 Ne and 40 Ti. 

The unit system h = c = 1 is used in the article. 



II. THREE-BODY RADIATIVE CAPTURE 

The derivations provided below are relatively trivial. 
Those of Section III Al can be found e.g. in Refs. [ll |23|. 
They are presented, however, in much details to provide 
unified notations, simplify the reading of the paper and 
to avoid any possible misinterpretation. 



A. Sequential capture 

The abundance Ya+i for the nucleus with mass number 
A + 2 due to the sequential two-proton capture reaction 
on the nucleus with mass number A is defined via three- 
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body reaction rate ( (see e.g. Ref. 

Y A+2 = (1/2) TVl p 2 (a ppn v) Y 2 Y A , (1) 

where p is density and iV^ is Avogadro number. The 
two-proton reaction rate is defined for sequential capture 
of protons (for example in 1]) by 



p. i 



This expression is a consequence of the rate equations for 
the resonance number i: 



yjji = N AP (a p , p v) i YpYA-T^YXl, 

Y A+ 2 =Y. N *P (0p,7«>< Y P Y M-l 



(i) 



(2) 



and the assumption about thermodynamic equilibrium 
for the intermediate resonant states V^i-i = 

The standard expression for cross section of the reso- 
nance reaction with entrance channel a and exit channel 
j3 is 

a( ' ~ k 2 2 (E - E R y + T 2 /4 (2Ji + 1)(2J 2 + 1) ' 

. (3) 

where J\ and J 2 are the total spins of incoming particles 
and J 2 r is the total spin of the resonance. 

In the case of intermediate capture into the narrow 
proton resonance number i (which also decays practically 
only via proton emission) T a — Tp 
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where Jj is the total spin of the initial (core) nucleus 
and J 2 R t i is the total spin of the resonance number i in 
the core+p system. The Boltzmann distribution over the 
relative motion momentum fc 12 = \/2mi 2 Ei 2 is 

w{k 12 ) = {2irm 12 kT)~ 3/2 exp[-E 12 /kT] , 

and we approximate the integral over the resonance pro- 
file as 



dE 



2n 



(e - Er) 2 + r 2 /4 r 



The reaction rate for the subsequent capture of the sec- 
ond proton on the system core+p in the resonant state 
number i and the following gamma emission is: 



A x + A 2 + A 3 



3/2 



2J F + 1 



2tt 



3/2 



(Ai + A 2 )A 3 ) 2{2J 2R . l + 1) \mkT y 

x e xp [-{E 3R -E 2R>i )/kT]r 1 r' P!i /r 3R , (5) 



where Jp is the total spin of the resonance E 3R in the 
corc+2p system. T' p i is the partial width for decay of 
this state into the binary channel (core+p) +p, where the 
core+p system is in the resonant state number i. It is 
easy to find out that the two-proton reaction rate for 
sequential proton capture (through narrow intermediate 
resonances) is 



1 pp~i 1 



A x + A 2 



A 3 \ 3/2 2J F 



AiA 2 A 3 J 2(2J/ + 1) \mkT J 

x exp [-E 3R /kT] (r 7 J2i r;,*) /r 3R . (6) 

The following features of the rate Eq. © should be noted: 

(i) The reaction rate does not depend on the number and 
properties of the intermediate states, but only on the sum 
of the proton widths for population of these states. 

(ii) In the most expected case of the sequential decay 
mode dominance for the three-body resonance E 3 r we 
have T 7 <C T 3 r, T 2p -C T 3R (T 2p is the width for the di- 
rect decay into 2p continuum, the process not proceeding 
via intermediate core+p resonances) and T 3 r = J2i Ip,i ■ 
So, the reaction rate depends only on gamma width of 
the three-body resonance. 

(iii) Eq. © shows that if there exist other significant de- 
cay channels for three-body resonance E 3 r, then T 3 r > 
J2 t T' p i and the production rate decreases. Such possible 
decay channel is, already mentioned, direct (not via reso- 
nances) decay of the three-body resonance E 3 r into two- 
proton continuum. Typically this process is suppressed, 
but there are cases where this process is not suppressed. 
It can be dominating, or even it can be the only possi- 
ble decay channel (no intermediate resonances). Such an 
opportunity is considered in the next Section. 



B. Direct capture 

The abundance Ya+2 due to the direct two-proton cap- 
ture reaction is defined via three-body reaction rate as 

Y A+2 = (1/2) N 2 A p 2 (a 2p!l v) Y 2 Y A . (7) 

To derive the cross section of the direct capture from the 
three-body continuum we use the S-matrix formalism for 
3^3 reaction. The plane wave for three particles can 
be decomposed over hyperspherical harmonics 1 1 j^ L '■ 

^ 3 W = exp[ikiri + ik 2 r 2 + ik 3 r 3 ] XS 1 m 1 Xs 2 m 2 Xs 3 m 3 
= exp[ik cm R cm + ikyY + ik x X\ XS 1 m 1 Xs 2 m 2 Xs 3 m 3 

(27T) 3 



exp[ zk cm R, cm ] 
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JM KLl m l v SS w 



E 



i J K +2{xp) 



9mm l m x m 2 m 3 > 
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where XSiMi are spin functions for the cluster number i. 
J JM denote the hyperspherical harmonic J LMl coupled 
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with spin functions of clusters to the total momentum J 

i^KUy ® x ssJjm , 
[\XSi ®XS 2 ]s x ®XSi]sm 3 ■ 



n-JM 

j kli~u.ss x 



XsS^Ms 



Variables k cm and R cm describe the center of mass mo- 
tion; {k^ky} and {X, Y} are conjugated sets of Jacobi 
variables for internal motion of the three-body system. 
Complete definitions of the hyperspherical variables and 
hyperspherical harmonics can be found e.g. in Ref. |2(1 |. 
The dependence on magnetic quantum numbers is 



JLSS 
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Following the same steps as in the two-body case we de- 
fine scattering amplitude and decompose it over hyper- 
spherical harmonics. Asymptotic form of the three-body 
WF (center of mass motion is omitted) is given by 



* 3 ( P - oo) = *r 

where 

fM 1 M 2 M 3 {Qp, fi*) 



ex P[^P] f (a o 

gy^ JM 1 M 2 M 3 (Hp, "> 
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X fMM L M 1 M 2 M 3 (^ l P^x) XsSxMs ■ 

3 scattering amplitude can be written as 



r J L S S x 
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(il p ,Sl x ) = exp[-wr/4] (2ir/x) 
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m 'l 

where 7 = {Ll x l y SS x }. In these equations angles fi^ 
point to the direction on the hypersphere, where the par- 
ticles have come from [the directions defined by momenta 
in the plane wave ijHJ]. The angles f2 p in the asymp- 
totic expression (JHJ define the vectors of particles and 
the energy distribution after collision. The cross section 
of 3 — > 3 scattering can be written as 



da(Q > 



SS X M S LM l 



n JM rJLSSx 



J LM L SM S J MM L Mi M 2 M 3 \"P 



and after integration over the angle fl p of the outgoing 
particles, summation over the projections of final total 
spin M and averaging over the projections of spins of the 
incoming clusters Mi 



<r(fi x ) = (2n/x) 5 G J SlS2S3 S S 



s ' 4: 6 



M L 



E 
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Coefficients Gg iS2 g are combinatorial factors 



C J 

u s 1 s 2 s 3 



2J + 1 



(2Si + l)(2S 2 + l)(2S 3 + l) 
The astrophysical production rate is given by 
x 



\CT2p,2 p V) 



a(fl x ) w(kik 2 k 3 ) d 3 kx d 3 k 2 d 3 k 3 . 



where the Boltzmann distribution for 3 particles is 

, , , eM-(Ei+E 2 + E 3 )/kT] 
W(klk2k3) = {2,mkTf/ 2 (A lA2 A 3 )^ ■ 

After transformation to hyperspherical variables 



d 3 fcid 3 /fc 2 d 3 fc 3 



A X A 2 A 3 



A 1 + A 2 + A 3 



3/2 



d 3 k cm dVlx x 5 dx , 



and after d 3 k cm and dfl^ integration (we should remind 
here that angles f2 x point to the directions of the incom- 
ing particles defined by vectors k,-) 



(0"2p,2pU)' 



{2nf 





" x 2 ' 


X 


j exp 


2mkT 


m 



7r(27rmfcT) 3 616 
x E K';'-S^y(S^'-S^')dx.(12) 

To get the inelastic part of the cross section for the suf- 
ficiently narrow resonance we should replace 0> US 



E « 7 '-s 

Ky,K J 7' 



K 'i'\Ux K 'i 



S 



T2 p r 7 



{E-E 3R f + Tl R /A 



(13) 



For the two-proton capture on the nucleus Ji to the final 
state Jf (assuming the small width of the three-body 
resonance) we obtain 



z / 2tt \ zj f + 1 
( CT 2 P ,7«> = ( Z^Ur ) 0^0 7_ I ^ ex P „ , 

The production rate for two-proton capture should be 
multiplied by squared density p 2 to provide the abun- 
dance, see Eq. {7}. Eq. lfTl|) is written in scaled Jacobi 
variables. The density in these variables can be expressed 
via density in the ordinary space as 



2J F 



1 



3-R 



kT 



T2pF 7 
T 



Ml 



2 _ 2 

/^scaled P 



Ai + A 2 + A 3 
A X A 2 A 3 



3/2 



The expression for production rate which can be used 
with expression for density in normal space (which is in- 
dicated by the absence of the prime symbol) is therefore 



(C2p/ 



i 1 + 1- : A 3 \ 3/2 2J F + 1 / 2tt \ ' 
2(2J/ + 1) \mkTj 

(15) 



A X A 2 A 3 

x exp [-E 3R /kT] r 2 pr 7 /r 3fl 
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Eq. i|15|) is absolutely the same as Eq. © except for the 
dependence on the decay width to the three-body contin- 
uum T 2p instead of decay widths to the resonant states in 
A+l system T' p t . We can draw the following conclusions 
here: 

(i) Eq. i|15|l is obtained in a very general assumptions 
about existence of the asymptotic © and analytical 
properties of the 3 — > 3 scattering S-matrix. We also use 
the fact that most of the states of interest are narrow. No 
other assumptions is made (e.g. in sequential formalism 
there is assumption about existence of the specific decay 
path). So, the direct capture (unlike sequential capture) 
is always possible. 

(ii) It is clear that Eqs. JfjJ) and i|15fl supplement each 
other and the total reaction rate is 



the asymptotic implied in derivations of Sec. Ill Al is 



(<7 2 p, 7 v) 



A 1 + A 2 



A 1 A 2 A 3 



3/2 



2 J F + 1 



2tt 
mkT 



exp 



E; 



:;./; 



kT 



r 



2p 



2(2J/ 
TT' 

t—ll p,t r 



1) 



7 • 



(16) 



(hi) In the most likely situation T 2p + T' p i 3> T 7 (and 
hence T 3 r = r 2p + Yli^'pi) tne total reaction rate de- 
pends only on gamma width of the three-body resonance. 
However, it is possible that the direct two-proton emis- 
sion is the only nuclear decay branch for the state. The 
width T 2p could be very small (smaller than the gamma 
width) in a relatively broad range of the three-body de- 
cay energies [2(J. In that case the reaction rate depends 
only on T 2p . 

(iv) Eq. I|16|) gives the formula for reaction rate, which is 
the same as one, known for a long time (see e.g. Eq. (20) 
in Ref. H^]), which was obtained by much easier means 
(namely, a complete thermal equilibrium and a detailed 
balance) than in this work. The result of our derivation 
here is clearer understanding of the fact that this formula 
already correctly and completely include both sequential 
capture and direct capture reactions. 

So, we see that for the resonant part of the reaction 
rate the sequential formalism treatment is overcompli- 
cated and incomplete. This is not a great issue in most 
cases, but there are situations, where it becomes impor- 
tant. The impact of the formalism on the rates of reac- 
tions of astrophysical interest is discussed in Sec. 11111 



C. Formal questions 

The derivations of the reaction rates above in Sections 
III Al and III Bl are quite schematic. They basically rely 
on assumptions about existence of definite asymptotics 
of the three-body problem, which could be not evident. 
They require if not a proof then at least some discussion. 

For sequential formalism we need that there exists a 
long-living resonance state in the X Jacobi subsystem 
(at energy E. x = k x /(2M x ) and with width T x ). Then 



* 3 ({X,Y} ► oo) = *f + —f(k x ) — f(k y ) 



Y 



For the direct capture the assumed asymptotic is 



* 3 (p ► oo) = *r + ^— f(n p , n x ) , 



(17) 



(18) 



where k 2 J(2M x ) + *J/(2M„) = x 2 /(2m) = E 3R . 

The expressions 11171 and (|18fl correspond to neutral 
particles, while we are speaking about nuclei Z — 8 — 20 
capturing protons. The typical densities for X-ray bursts 
and processes in novae are 10 3 — 10 6 g/ccm [28|. For such 
densities the average distances between protons are about 
6 x 10 2 — 2 x 10 4 fm and for characteristic temperatures 
the Debye screening radii are 3 x 10 3 — 5 x 10 5 fm. Beyond 
these radii we have a formal right to use the forms (|17|l 
and JTHJl. 

It should be understood that in a very formal sense the 
asymptotic <|171 is not valid. In the limit of infinite dis- 
tance a long-living two-body state finally decays and the 
asymptotic (|T7|) should be replaced by (JTSJ) . However, 
from practical side the separation of asymptotics (|17|) 
and (|18fl is reasonable. A nice example (also relevant to 
the further discussion) of the coexistence of the three- 
body and binary asymptotics is the decay of 9 Be 5/2~ 
state at 2.429 MeV. The branchings to the three-body 
channel Br(3) and binary channel Br(2) are compara- 
ble (Br(3) ~ 0.93 - 0.95 and Br{2) ~ 0.07 - 0.05 [U), 
which is experimentally well observed (see, for example, 
Ref. H3]). In the case of binary decay the average flight 
distance of the 8 Be g.s. resonance (r^ = 6.8 eV) is around 
10 6 fm. Thus it is clear that for some practical purposes 
the assumption of Eq. (|17fl is justified. It is also clear 
that the broader is the resonance in the X subsystem, 
the faster a transition from Eq. (|17fl to Eq. I|18|) happens. 
For example, for the width of the intermediate resonance 
T x around 100 keV (and typical E 3 r — 1 MeV) the av- 
erage flight distance of this resonance is around 100 fm. 
This is much smaller than the typical distance between 
protons in the stellar media and then usage of Eq. (|17|) 
and chemical balance description by the Eqs. loose 
sense. 

The asymptotics (|17f) and l|181 are located in the same 
space and could have the same quantum numbers. How- 
ever, we can can speak about orthogonality of these 
asymptotics in definite sense and then treat currents as- 
sociated with them independent. Only this assumption 
makes possible the separate treatment of sequential and 
direct decay channels in Sections III Al and III Bl The 
asymptotic (I17f) is typically localized in a very small part 
of the phase space. Formally this corresponds to the fact 
that the hypersphcrical scries for binary channel (at given 
hyperradius) is very long (with many significant terms), 
while for asymptotic (|18fl we can expect that only the 
lowest hyperspherical harmonics in the decomposition are 
significant. Thus the sequential and direct channels are 
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practically orthogonal on the hypersphere of a large ra- 
dius. With increase of the radius this "orthogonality" 
(called asymptotic orthogonality) becomes better, until 
the effect of the X subsystem decay becomes important. 
Again, on the example of the 5/2 _ state of 9 Be, the level 
of overlap between the direct and sequential decay chan- 
nels in the momentum space can be estimated as 



Br(3) Br{2) 



8 B< 



E 3 r(5/2-) 



0.06 



6.8 eV 
764.1 keV 



5 x 10" 



which is clearly a very small value. When there is no such 
a reliable separation between channels i|17fl and (|18fl any 
more (the intermediate resonances are too broad) we have 
a formal right to speak only about asymptotic (|18|) . 

To finalize this discussion, the derivations in Sections 
III Al and III Bl were done as if only one type of the asymp- 
totic exists. From formal point of view only asymp- 
totic of Eq. (|18|l exists. For practical purposes either 
(i) only asymptotic of Eq. I|18fl exists or (ii) both asymp- 
totics Eqs. i|17fl and l|18fl are present simultaneously. The 
asymptotic (|18(l exists in the three-body problem uncon- 
ditionally, while the existence of (|17fl is subject to avail- 
ability of sufficiently narrow intermediate resonances for 
the decay. In the case (ii), the regions of the domina- 
tion of each asymptotic are separated by a complicated 
surface in the phase space. Some discussion of the rel- 
evant questions can be found, for example, in Ref. |25j |. 
So, the phase space integration in Sections III Al and III Bl 
should have been done not over all the space, but over 
regions of validity for each type of asymptotic. This is 
clear, however, that this imperfection does not influence 
the final result. The reasons are that (i) the contribution 
of the asymptotic of the selected kind in the phase space 
outside the region of its domination is typically negligible 
and (ii) anyhow we are interested in the contribution of 
the both kinds of asymptotics simultaneously. 




T (GK) 

FIG. 1: Reaction rate for lj O(2p,7) 17 Ne reaction. Solid curve 
shows calculations of this paper. Gray curves indicate bound- 
aries due to uncertainties in the input (e.g. the upper gray 
curve is obtained with parameter set r™ ax from Table [HJ . 
Dashed and dotted curves show full result from Ref. £jj and 
resonance contribution to it respectively. 



TABLE I: Reaction rates multiplied by N\ (in cm 6 /s) for 
15 0(2p,7) 17 Ne reaction. 



T (GK) 


Ref. [1] 


This work 


This work upper 


0.3 


2.9 x 


1Q -23 


4.9 x 


1Q -19 


2.9 x 


1Q -14 


0.5 


6.0 x 


10 -is 


2.1 x 


10 -15 


1.3 x 


IO" 12 


0.6 


1.2 x 


10 -16 


2.8 x 


1Q -14 


2.8 x 


IO" 12 


0.8 


5.6 x 


10 -15 


6.3 x 


10 -13 


6.9 x 


1Q -12 


1.0 


5.0 x 


1Q -14 


3.5 x 


lO" 12 


1.3 x 


io- 11 


1.5 


1.1 X 


1Q -12 


2.5 x 


HP 11 


3.8 x 


10 -n 


2.0 


6.0 x 


1Q -12 


5.1 x 


io- 11 


6.8 x 


io- 11 


3.0 


4.3 x 


10" 11 


7.6 x 


io- 11 


1.1 X 


10 -io 


5.0 


2.5 x 


10 -io 


7.3 x 


1Q -ii 


2.3 x 


10 -io 



III. DISCUSSION 



A. 15 0(2p,7) 17 Ne reaction 

The results of rate calculations for this reaction are 
shown in Fig. ^ and in Table [I] They differ significantly 
from the results of Ref. [l| (shown in Fig. ^ by dashed 
and dotted curves). For the temperature range of astro- 
physical interest (~ 0.3 — 3 GK, see [2^, for example) 
the expected increase of the rate, compared to Ref. [l| , is 
up to 4 orders of the magnitude, while maximal possible 
increase is up to 9 orders of the magnitude. 

The reasons of difference are evident from Table [H] 

(i) The level scheme of 17 Ne has been somewhat updated 
(see e.g. Ref. [2^) since the the work pj had been written. 

(ii) The use of Eq. ifTB^I includes the first 3/2 _ excited 
state of 17 Ne into treatment (it was omitted in Ref. @], as 
there is no sequential capture path to this state) . The im- 
portant difference of the situation with this state from the 



others is that gamma width of this state is known to be 
much larger than the 2p width and the reaction rate Eq. 
(|16fl is entirely defined by the 2p width for the simultane- 
ous two-proton emission. At the moment there exist two 
theoretical calculations of this width T 2p = 4.1 x 10~ 16 
MeV [13, T 2p = 3.6 x 10~ 12 MeV [2l|, and a quite re- 
laxed experimental lower lifetime limit of r > 26 ps [l6| 
(which corresponds to the width T 2p < 2.5 x 10 -11 MeV). 
Values T 2p from Refs. 0] and are used to estimate, 
respectively, the lower and the upper boundaries for the 
band of expected values of the rate (see Fig.^l. The res- 
onance contribution of this state is dominating the rate 
in the temperature range 0.05 — 0.35 GK if we take the- 
oretical 2p width from Ref. ^?|], an d up to 1.2 GK if we 
consider the experimental limit. 

(iii) In paper pj the gamma widths for 17 Ne were taken 
from studied transitions in the isobaric mirror partner 
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TABLE II: Resonance parameters of 1 Ne states used in the 
15 0(2p,7) 17 Ne reaction calculations. The column T™ shows 
the upper experimental limit for T 7 values, and F™ m the val- 
ues used for lower limit estimates. The listed set of states is 
sufficient for rate calculations up to 5 GK. The widths and 
branchings, which are not specially discussed are from |2q. 



State 


Ref. [U 




This work 




r 


E (keV) 


T 7 (eV) 


E (keV) 


T™ in (eV) 


r^ nax (eV) 


3/2- 






1288 


"4.1 x 10 -10 


b 2.5 x 10 -5 


5/2" 


1907 


6.0 x 10- 5 


1764 


c 1.7 x 10- 3 


c 2.0 x 10- 3 


1/2+ 


1850 


1.6 x 10~ s 


1908 


1.6 x 10- 5 


f 2.1 x 10" 4 


5/2+ 


2526 


2.0 x 10 _s 


2651 


d 9.0 x 10" 6 


f 9.9 x 10" 5 


3/2- 


3204 


0.022 


3204 


e 0.019 


f 0.19 



"This is a 2p width (see Eq. 1161 : the gamma width is dominating 
the decay of this state). This value is calculated theoretically in 
Ref. [JJ. 

This is a 2p width. This experimental limit on 2p width is found 
in Ref. |Jj|. 

This value is calculated from B(E2) = 124(18) e 2 fm 4 given in 
Ref. E3. 

d This value is partial width (45% branching) into the l/2~ ground 
and the first excited 3/2- states of 17 Ne. The gamma transition 
to 1/2+ state returns the system into 2p continuum. 

e This value is partial width (88% branching) into the ground and 
first excited states of 17 Ne. 

^These values are assumed (in analogy with more than order of 
magnitude increase for 5/2" state from column 3 to column 5). 



TABLE III: Resonance parameters of Mg states used in the 
18 Ne(2p,7) 20 Mg reaction calculations. 



State 


This work lower 


This work upper 


r 


E (keV) 


T 7 (eV) 


E (keV) 


T 7 (eV) 


4 + 


3570 


2.1 x 10" 4 


3451 


5.6 x 10 -4 


2+ 


4072 


1.3 x 10 -3 


3857 


8.9 x 10" 2 


o + 


4456 


1.4 x 10 -3 


4317 


1.4 x 10 -3 


4 


4850 


2.6 x 10- 3 


4699 


2.6 x 10" 3 


2 3 + 


5234 


2.9 x 10 _1 


4978 


2.9 x 1Q- 1 



TABLE IV: Reaction rates multiplied by N\ (in cm 6 /s) for 
18 Ne(2p,7) 20 Mg reaction. 



T (GK) Ref. [1] This work lower This work upper 



0.3 


4.4 x 10" 


■21 


9.3 


X 


10" 


28 


2.5 x 


10" 


■28 


0.5 


3.3 x 10" 


-17 


4.3 


X 


10" 


-20 


1.8 x 


10" 


■18 


0.6 


4.8 x 10" 


-16 


3.0 


X 


10" 


■18 


8.3 x 


10" 


■17 


0.8 


1.8 x 10" 


-14 


5.1 


X 


10" 


16 


9.5 x 


10" 


■15 


1.0 


2.4 x 10" 


-13 


9.6 


X 


10" 


■15 


1.8 x 


10" 


■13 


1.5 


6.8 x 10" 


-12 


3.7 


X 


10" 


13 


1.1 X 


10" 


■11 


2.0 


9.7 x 10" 


-11 


2.0 


X 


10" 


■12 


8.0 x 


10" 


■11 


3.0 


4.8 x 10" 


-10 


1.5 


X 


10" 


■11 


4.4 x 


10" 


■10 


5.0 


2.4 x 10" 


-9 


1.4 


X 


10" 


■10 


1.2 x 


10" 


-9 



17 N. Recently the decay of the first excited states of 
17 Ne (3/2~, 5/2-) has been studied via intermediate en- 
ergy Coulomb excitation of a radioactive 17 Ne beam on 
a Au target In this paper the transition matrix 
elements B(E2,l/2~ -> 3/2") and B(E2,l/2 - -> 5/2") 
have been deduced. We use the deduced B(E2) value 
from Ref. to calculate the gamma width of the 5/2~ 
state. The result is shown in Table ITU This width ap- 
pears to be about 30 times larger that the corresponding 
width of the mirror state in 17 N. This is, probably, con- 
nected to the fact that at the proton-rich side the number 
of the protons contributing gamma transitions is larger 
and these protons are situated at larger distances com- 
pared with tightly bound protons in 17 N. This situation 
is also expected for the other states in 17 Ne (compared 
to the states in 17 N), which is reflected by an order of 
the magnitude increase of the other widths for estimates 
of the upper limits (column IT ax in Table |TTJ. 



B. 18 Ne(2p,7) 20 Mg reaction 

For the 18 Ne(2p,7) 20 Mg reaction there are no three- 
body states which were not taken into account in Ref. 
1], so, no significant update of the rate is expected here. 
However, the level scheme and gamma widths are not 
known experimentally for this nucleus and this should 
be reflected in the rate calculations 



In the cases of 2p capture into 17 Ne and 40 Ti the 
gamma widths from mirror isobaric partners were used 
in Ref. [jj]. In contrast, for capture into 20 Mg the system- 
atics values were utilized. The theoretical B(E2) values 
for some low-lyi ng s tates in 20 Mg have been calculated 
recently in Ref. jajj. The £(E2,4+ -> 2+) was found 
to be 28.2 or 11.6 e 2 fm 4 (for V2 and MN forces respec- 
tively) and £(£2,2^ -> 4 ") was found to be 2.9 or 1.9 
e 2 fm 4 . These reduced probabilities give gamma widths 
5.6 x lfj- 4 or 2.3 x 10~ 4 eV for 4+ state (which is compa- 
rable to value 2.1 x KT 4 eV used in Q) and 2.0 x lO" 3 or 
1.3 x lO -3 eV for 2 J state (which is significantly less than 
8.9 x 10" 2 eV used in 0). We combine the larg est and 
the lowest gamma widths from Refs. and |30j] to esti- 
mate the upper and the lower boundaries for the rate (see 
Table ITTT)l . To incorporate in this estimate the sensitivity 
to the level scheme we also use for the lower estimate the 
energies of the states from 20 O. The distance between lev- 
els here is expected to be somewhat larger than in 20 Mg 
and the reaction rate thus should further decrease. 

The results of calculations are shown in Table llVl The 
upper boundary in our calculations is in a good agree- 
ment with results of Q (factor of two) at T > 0.8 GK. It 
was shown in that below 0.8 GK the nonresonant con- 
tribution to the reaction rate dominates, which explains 
the discrepancy in Table HVl at low temperatures. 
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TABLE V: Resonance parameters of Ti states used in the 
38 Ca(2p,7) 40 Ti reaction calculations. In the "Type" column 
the type of width is specified, which defines the contribution 
of the state to reaction rate. The "Lower" set of widths is 
used in the calculations with S2 P = 1.370 MeV, while the 
"Upper" set with S 2p = 1.582 MeV. 



,F E (keV) Type "Lower" L (eV) "Upper" T (eV) 



o 2 + 


2121 


T2p 


io- 12 


io- 15 




■4 


2524 


r P 


1.0 x 10~ 5 


1.0 x 10" 


5 


*t 


2892 


r 7 


2.0 x 10 -4 


2.0 x 10" 


3 


2.1 


3208 


r 7 


1.0 x 10~ 2 


1.0 x 10" 


1 



TABLE VI: Reaction rates multiplied by N\ (in cm 6 /s) for 
38 Ca(2p, 7 ) 4() Ti reaction. 



T (GK) Ref. [1] This work lower This work upper 



0.3 


2.1 


X 


10 


25 


7.0 


X 


10 


28 


2.3 


X 


10 


24 


0.5 


1.0 


X 


10" 


19 


7.8 


X 


10" 


21 


1.1 


X 


10" 


18 


0.6 


2.4 


x 


10" 


18 


4.0 


X 


10" 


19 


3.1 


X 


10" 


17 


0.8 


1.1 


x 


10" 


16 


5.3 


X 


10" 


17 


3.0 


X 


10" 


15 


1.0 


1.3 


x 


10" 


15 


1.2 


X 


10" 


15 


7.1 


X 


10" 


14 


1.5 


7.0 


x 


10" 


14 


1.3 


X 


10" 


13 


6.0 


X 


10" 


12 


2.0 


5.2 


x 


io- 


13 


1.5 


X 


10" 


12 


5.1 


X 


io- 


11 


3.0 


3.0 


x 


10" 


12 


1.4 


X 


10" 


11 


3.2 


X 


10" 


10 


5.0 


8.0 


x 


10" 


12 


5.0 


X 


10" 


11 


8.2 


X 


10~ 


10 



C. 38 Ca(2p,7) 40 Ti reaction 

For the 38 Ca(2p,7) 40 Ti reaction one three-body state 
was omitted in Ref. According to the isobaric sym- 
metry there should be a 0^" state located at about 2.121 
MeV excitation energy. The two-proton separation en- 
ergy used in Ref. is S 2p = 1.582 MeV. Another esti- 
mate (e.g. [3l]]) is S2 P = 1.370 MeV. In the first case the 
2p emission energy for 0^" state is 539 keV and (following 
Refs. 0,|2(j) the two-proton width can be estimated as 
about 10~ 21 MeV. In the second case the 2p energy is 751 
keV and the estimated two-proton width is around IO -18 
MeV. For other states we use parameters from Ref. pj 
(see Table Ivjl . which mainly come from the isobaric mir- 
ror partner 40 Ar. To estimate the upper boundary for 
the reaction rate we increase the gamma widths of 4^ 
and 2^ states by an order of the magnitude. As we have 
already discussed, one could expect a significant increase 
of the gamma widths when we come to the proton-rich 
mirror partner. To estimate the sensitivity to the level 
scheme (which is not known for 40 Ti) we use the smaller 
2p separation energy S 2 p = 1-370 MeV for the estimate of 
the lower boundary and S2p — 1.582 MeV for the upper 
boundary. Again, as in the case of 20 Mg, the increase of 
the state energy above 2p threshold leads to decrease of 
the corresponding reaction rate. For that reason we use 
the larger 2p width for estimate of the lower boundary 
(Table Ivl line 1). The larger two proton width corre- 
sponds to the case of larger energy of states above 2p 
separation threshold. 

The results of calculations for 38 Ca(2p,7) 40 Ti are given 
in Fig. [3 and Table IVTl Our results are somewhat larger 
(1-2 orders of magnitude) than results of [l| for temper- 
atures T > 1 GK. They more or less overlaps at lower 
temperatures. The effect of inclusion of Oj" 2.121 MeV 
state can be seen in Fig. the range between upper and 
lower boundaries shrinks at T < 0.35 GK. This happens 
because the contribution of the 0^" state is much larger 
in the "Lower" parameter set, which otherwise provides 
a smaller reaction rate. 




0.1 1 5 

T (GK) 



FIG. 2: Reaction rate for 38 Ca(2p,7) 40 Ti reaction. Solid curve 
shows the result from |l|. Gray curves indicate upper and 
lower boundaries for our results (see Table Wl\ due to existing 
uncertainties in the input. 



D. 4 He(na,7) 9 Be reaction 

The stellar reaction rate for 4 He(na, 7 ) 9 Be process has 
been studied several times in the recent years 0, 0, IE S 
The results are in overall agreement, except for the 
latest paper 0- In this work the rate is obtained which 
is significantly higher (for temperatures T > 3 GK) than 
the rates in the previous studies. 

In our studies here we have found that the sequen- 
tial formalism underestimate the reaction rate only if the 
width of the state for direct decay into continuum is dom- 
inating (see Sectioning). The low- lying (E < 3 MeV) 
9 Be states typically have strong 8 Be+n decay branch- 
ings. Only the 5/2~ 2.429 MeV state is an exception: 
the branching to the three-body channel is 93 — 95% pfj) . 
The gamma width of this state is 0.091 eV 0. The re- 
sults of our calculations are shown in Fig. [3] and Table 
IVIII In these calculations we use a version of Eq. (|Pd|) 
without an assumption about narrow widths of the res- 
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TABLE VII: Reaction rates multiplied by N% (in cm 6 /s) for 
4 He(na,7) 9 Be reaction. 



FIG. 3: Reaction rate for 4 He(na,7) 9 Be reaction. The solid 
curve shows our result (see also Table IVTTl . The dashed curve 
is the same, but without 5/2~ 2.429 MeV state contribution 
(this coincides with result of Ref. 0). The dash-dotted curve 
is contribution from the near threshold l/2 + state. The dot- 
ted curve shows calculations Ref. 01 . 



onanccs and the capture cross section is parametrized as 
in Ref. (with exception that 5/2 ~ state is included). 
The results obtained are in a very good agreement with 
0. The increase of the rate due to addition of the 5/2 ~ 
state is 11% at most in the temperature range up to 10 
GK. This small change is connected with comparatively 
small gamma width of this state: the gamma widths of 
the other states in the capture cross section parametriza- 
tion used in || are around 0.45 — 0.9 eV. So, the un- 
certainty of the reaction rate due to uncertainties of the 
experimental data found in |]| is significantly larger than 
correction connected with 5/2 ~ state (see Table fyTTl . 

The mentioned experimental uncertainty could be even 
larger than it was inferred in Ref. Q. The analysis, 
provided in Ref. in the framework of the semimicro- 
scopic model, demonstrated that the older photodisintc- 
gration data for 9 Be (32,113,13 cou ld be more preferable 
than the more up-to-date results [35j (on which, e.g. the 
parametrization of the cross section used in Ref. |8( is 
based). The reaction rate found in paper (as well as 
in the early work Q) is around 35% larger than the rate 
in Ref. @. 

Paper [j| is generally dedicated to the R-matrix anal- 
ysis of the /3-delayed particle decay of 9 C via the excited 
states in 9 B. The authors utilize the R-matrix parame- 
ters obtained in the decay studies of 9 B for the caption 
calculations in 9 Be. The reaction rate calculated in this 
work is consistent with the other results at low temper- 
atures, but is qualitatively different at T > 3 GK (see 
Fig. |3 dotted curve). The rise of the reaction rate at 
higher temperatures is connected, according to Q, with 
contribution of sequential capture of a-particle on the 
broad ground state of 5 He. Such capture path has never 
been considered elsewhere. It should be noted that in the 



T (GK) 



Ref. 



Ref. [8] upper This work 



2 
4 
6 
8 

10 



1.80 x 10" 7 
5.48 x 10~ 8 
2.88 x 10 -8 

1.81 x 10~ 8 
1.23 x 10~ 8 



2.20 x 10 -7 
6.99 x 10~ 8 
3.83 x 10" 8 
2.47 x 10" 8 
1.70 x 10~ 8 



1.83 x 10~ 7 
5.94 x 10" 8 
3.20 x 10~ 8 
2.02 x 10~ 8 
1.36 x 10~ 8 



framework of sequential formalism this is a valid ques- 
tion: how narrow should be the intermediate state, to be 
considered within this formalism. Really, the sequential 
formalism is evidently correct in the limit of infinitely 
narrow intermediate state. However, in the other limit 
(an infinitely broad state), we have just nonresonant con- 
tinuum and the sequential formalism should fail at some 
point. This issue is qualitatively discussed in Section llI CI 
Our work resolve this question in a very natural way: we 
state that contributions of different sequential and three- 
body channels should add up in a way, which makes their 
relative contributions unimportant. So, inclusion of cap- 
ture via 5 He into formalism should not lead to any signif- 
icant changes (compared to conventional sequential cap- 
ture via 8 Be g.s.), until there exist states with dominating 
three-body decay branch (which is not accounted in se- 
quential formalism) and large gamma widths. No such 
states are known in the energy range of interest. The 
reaction rate from Ref. can be reproduced within our 
formalism only if we assume the gamma width for the 3 
MeV state in 9 Be to be about 15 eV and also assume one 
more state at about 5 MeV with gamma width above 1 
keV. Such assumptions are quite unrealistic. 

Unfortunately, there is an evidence for problems in 
Ref. 0, which probably have leaded to the discussed 
strange result. In Eq. (32) of this work the penetrability 
is present in the first power, while it should be in the sec- 
ond (as we speak about elastic cross section). Possibly 
this is the reason of the qualitatively incorrect behaviour 
of the intermediate population values (see Fig. 8 in Ref. 
Q). For example, the population (a(E)v/T(E)) for 5 He 
g.s. should decrease as T at low temperature. In Fig. 8 
of Ref. this value has a rapid rise at low temperature. 
Using Eq. (32) from Ref. "as is" one gets behaviour 
T -1 / 2 at low T in agreement with this figure. 

So, the difference of our approach found in the case 
of 4 He(na,7) 9 Be reaction is not significant. It is much 
smaller than the other uncertainties (see Refs. and 
8J). However, our formalism excludes such a possibility 
as an importance of broad intermediate 5 He state |9j. 
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IV. CONCLUSION 

We use the formalism based on the S-matrix for 3^3 
scattering to derive the reaction rate for the three-body 
resonant radiative capture. This derivation makes espe- 
cially evident that (i) all the three-body states should be 
included in the treatment (even if there is no opportu- 
nity of a sequential capture to the state), (ii) the detailed 
knowledge of the intermediate states is unnecessary to 
calculate the resonant rates and (iii) only the knowledge 
of particle and gamma widths for the three-body states 
is needed to calculate the resonant rates (not the relative 
contribution of direct and sequential mechanisms). 

This formalism, together with the modern results on 
2p and 7 widths of 17 Ne states, allows us to update signif- 
icantly the capture rate for the 15 0(2p,7) 17 Ne reaction. 
The updated rate is up to 4-9 orders of the magnitude 
larger (in the temperature range of astrophysical inter- 
est). The experimental derivation of the 2p width of the 
first excited state in 17 Ne is found to be very important 
for refining this rate. The 38 Ca(2p,7) 40 Ti reaction rate 
has also got a considerable increase. Thus the conclu- 



sions about importance of the 2p capture reactions could 
possibly be more optimistic than in Ref. [| . We also dis- 
cuss the impact of our approach on the 18 Ne(2p,7) 20 Mg, 
and 4 Hc(na,7) 9 Be reaction rates. Our studies emphasize 
the importance of better gamma width information for 
2p capture rates (experimental or theoretical, if the first 
is not available). 

The studies of this work are restricted to resonant re- 
actions (and correspondingly to relatively high tempera- 
tures). We are planning to perform accurate three-body 
studies of the nonresonant contributions in the forthcom- 
ing paper. 
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